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ABSTRACT 

In  order  to  simplify  the  construction  of  the 
computational  grid  about  complex  geometries  and  to 
enable  the  numerical  simulation  of  relative  motion 
between  components,  an  overset  grid  assembler 
SUGGAR  has  been  integrated  with  an  existing 
RANS  flow  solver  UNCLE.  This  modified  flow 
solver  with  overset  capability  allows  different  flow 
regions  to  have  their  own  separate  grid  that  is 
adaptive  to  the  particular  feature  of  that  region.  The 
complete  flow  field  is  then  made  up  of  a group  of 
overlapping  sub-grids.  A grid  assembler  is  then  used 
to  provide  the  flow  solver  with  the  connectivity 
information  among  the  sub-grids  and  the 
interpolation  stencils  of  the  boundary  nodal  points. 
To  demonstrate  the  flexibility,  feasibility  and 
accuracy  of  the  overset  grid  scheme,  numerical 
simulations  of  flow  over  ONR-Body-1  were 
performed. 

ONR-Body- 1 is  a frilly  appended  submarine 
model.  It  was  built  with  the  intention  to  study  the 
complicated  flow  field  about  a submarine  during 
maneuvering.  It  was  tested  in  the  tow  tank  at 
NSWCCD.  Three-dimensional  velocity  components 
around  the  body  were  measured  with  the  SPIV 
technique,  the  force  and  moment  coefficients  were 
obtained  through  block  gages. 

Comparisons  of  simulation  and  experiment 
and  discussion  are  presented  in  the  report. 

NOMENCLATURE 

Aref  Reference  area,  (Lref ) 2 

„ P 

Up  pressure  coefficient  — 

KpVr^f 


F 

Force  on  body,  in  x,y  or  z direction 

k',m’,n' 

Moment  coefficients 

M 

X pYref  Lref  Aref 

roll,  pitch  and  yaw. 

M 

Moment  about  the  center  of  gravity, 
roll,  pitch  or  yaw 

Lref 

Reference  length,  hull  length 

P 

Static  pressure 

Re 

Reynolds  number  (VrefLref  )/t) 

vref 

Reference  velocity  , ship  speed 

x,y,z 

Cartesian  coordinate  system, 
origin  located  at  the  forward 
perpendicular  along  the  center  line  of 
the  hull 

X , Y ,Z 

F 

Force  coefficients  , 

Vl  PYrgf  Aref 

in  x,y,  z directions 

y+ 

Turbulence  coordinate,  d-Jxw/p/t) 

u 

shear  stress  at  wall 

p 

Density 

t) 

kinematic  viscosity 

d 


distance  from  the  wall 


INTRODUCTION 

Finite  difference  and  finite  volume  methods  are  well 
developed  for  the  approximate  solutions  of  a system 
of  partial  differential  equations,  such  as  Reynolds- 
averaged  Navier-Stokes  equations.  Commonly,  the 
solution  domain  is  first  discretized  with  inter- 
connected grid  points  that  are  well  ordered.  The 
structured  nature  of  the  grid  points  makes  it  easier  to 
implement  the  high-order  approximation  of  the 
differential  equations.  In  the  case  of  solving  viscous 
flow  problem,  near  the  solid  surface,  the  grid  can  be 
easily  clustered  so  that  the  boundary  layer  profile  can 
be  sufficiently  resolved.  However,  for  a more 
complex  geometry,  where  a complicated  flow  pattern 
is  anticipated,  a highly  skewed  grid  is  difficult  to 
avoid.  In  general,  the  solutions  supported  by  skewed 
grids  are  less  accurate.  Generating  an  unstructured 
grid  system  for  complex  geometries  is  relatively 
easier  but  not  necessarily  trivial.  Besides,  because  of 
the  unstructured  grid  system  is  not  well  ordered, 
higher-order  approximation  is  not  easy  to  enforce 
[Steger,  Benek,  1987].  An  alternative  for  dealing 
with  the  complex  geometries  is  the  use  of  a 
composite  grid  system,  in  which  individual  structured 
grids  representing  different  regions  are  overset  on 
each  other.  Such  a grid  system  also  becomes  very 
desirable  when  treating  the  components  that  are 
moving  with  respect  to  each  other.  The  connectivity 
infonnation  among  the  grids  and  the  interpolation 
stencils  for  the  interface  boundary  nodes  can  be 
identified  and  determined  for  the  use  of  the  flow 
solver.  Clearly,  the  overset  grid  approach  requires 
additional  effort  for  grid  assembling,  but  in  return, 
the  process  of  construction  of  the  grid  becomes 
simplified,  especially  when  complex  geometries  are 
involved.  In  addition,  there  is  more  control  of  grid 
quality. 

Motivated  by  those  potential  benefits, 
Steger,  Dougherty  and  Benck  [Steger,  Dougherty, 
Benck,  1983]  first  introduced  the  overset  grid 
concept  in  1983.  Since  then,  much  research  efforts, 
especially  those  provided  by  the  Integrated  Wing 
Design  (IWD)  element  of  the  NASA’s  Advanced 
Subsonic  Technology  (AST)  program,  were  directed 
to  the  development  of  the  overset  CFD  software. 
More  recently,  a successful  compressible  flow 
simulation  of  a Boeing  777  aircraft  in  landing 
configuration  was  reported  [Rogers,  Roth.,  Cao, 
Slotnick,  Whitlock,  Nash,  Baker,  2000],  For  this 
particular  computation,  22.4  million  grid  points 
within  79  overset  zones  were  used.  To  support  and 
facilitate  the  computations  of  such  high  degree  of 
complicity,  several  grid-assembling  utilities  have 
been  developed:  (1)  OVERGRID,  to  generate  and 
manipulate  the  overset  grids  [Chan,  Gomez,  Rogers, 


Buning,  2002.],  [Chan,2002],  (2)  PEGASUS,  to 
perform  hole-cutting  and  to  provide  the  connectivity 
information  among  overlapping  grids  [Suhs,  Rogers, 
Dietz,  2002],  (3)  FOMOCO,  to  integrate  force  and 
moment  on  the  solid  surfaces  that  are  covered  with 
overlapping  grids  [Chan,  Buning,  1996],  Those  utility 
programs  have  become  workhorses  for  many  CFD 
applications  in  NASA,  universities  and  industries. 

In  the  field  of  incompressible  flow,  Kim 
simulated  rotor  blade  tip  vortex  cavitation  with  a 
overset  grid  approach  [Kim,  2002],  A dense  and  high 
quality  grid  is  placed  around  the  blade  tip  region  to 
resolve  the  fine  flow  structure  of  the  tip  vortex. 
Wilson  et  al  used  the  same  method  to  simulate  the 
large  amplitude  motion  and  maneuvering  for  surface 
ship  [Wilson,  Carrica,  Hyman,  Stem,  2004],  There, 
the  ability  to  use  the  local  refinement  for  the  Kelvin 
wave  system  with  the  overset  grid  is  illustrated.  The 
details  of  the  method  are  described  in  a reference  by 
Peterson  et  al  [Paterson,  Wilson,  Stem,  2004],  In 
essence,  this  incompressible  RANS  solver 
CFDSHIP-IOWA  is  based  on  pressure  implicit  with 
split  operator  (PISO)  algorithm.  The  grid  assembling 
was  accomplished  with  the  utility  programs 
developed  at  NASA  as  described  earlier. 

It  is  reasonable  to  conclude  that  the  Chimera 
grid  embedding  technique  is  desirable  and  mature  for 
CFD  applications.  The  fact  that  the  grid  assembling 
process  can  be  modularized,  makes  adding  the 
overset  capability  to  a flow  solver  extremely  feasible. 
It  is  a logical  and  natural  way  to  bring  the  existing 
and  validated  flow  solver  to  a new  and  higher  level  of 
application.  R.  Noack  has  developed  an  interface 
routine  DiRTlib  that  makes  the  process  of  adding  the 
overset  capability  to  a flow  solver  extremely 
transparent  [Roack,  2005],  The  DiRTlib  is  installed 
in  the  flow  solver  UNCLE,  this  new  version  with 
overset  capability  is  termed  as  UNCLE-overset.  A 
similar  endeavor  has  been  carried  out  at  PSU,  and 
some  numerical  results  were  reported  [Boger,  Dryer, 
2006],  This  report  describes  and  discusses  the  flow 
solver,  the  grid  assembler  and  the  numerical 
simulations  performed  for  the  validation  purpose. 
New  software  USURP  designed  by  Boger  to  integrate 
the  forces  and  moments  on  overset  grids  is  presented 
in  the  report. 

The  FLOW  SOLVER  - UNCLE-overset 

UNCLE  was  developed  at  Mississippi  State 
University  for  solving  the  three-dimensional 
incompressible  Reynolds-average  Navier-Stokes 
equations  [Taylor,  1991],  Artificial  compressibility  is 
used  to  couple  the  pressure  and  velocities  [Chroin, 
1967],  with  upwind  wind  flux-difference  splitting 
approximation  for  the  convective  terms  [Roe,  1981] 


and  second  order  central  differencing  approximation 
for  the  diffusive  terms.  With  this  approach,  die 
RANS  equations  are  treated  like  a system  of 
hyperbolic  equations.  UNCLE  has  accumulated 
certain  amount  validation  throughout  the  years.  To 
enhance  its  accuracy,  a fifth-order  stencil  was  added 
to  approximate  the  convective  terms  [Yang,  Jiang, 
Chesnakas,  Jessup,  2004],  To  perform  the  overset 
grid  calculation,  some  modifications  about  the  matrix 
of  the  system  of  equations  of  the  solver  are  needed 
for  the  fringe  and  the  hole  points:  (1)  set  the  off- 
diagonal  elements  to  zero,  (2)  set  the  diagonal 
elements  to  unity,  and  (3)  set  the  residuals  to  zero.  As 
a result,  the  solutions  on  those  points  are  not  updated 
during  iterations.  Interface  subroutines  are  also 
implemented  to  UNCLE  to  facilitate  the  overset 
calculation. 

For  turbulence  modeling,  there  are  two 
options  in  the  code:  the  k-£  and  q-co  models.  The 
results  presented  here  are  based  on  the  latter  model. 

THE  GRID  ASSEMBLER  - SUGGAR 

The  overset  grid  utilizes  a set  of  overlapping  grids  to 
discretize  the  domain  and  uses  interpolation  at 
appropriate  points  to  couple  the  solution  on  the 
different  grids.  Any  points  that  lie  outside  the  domain 
of  interest,  for  example,  the  grids  discretizing  the 
submarine  hull  that  lie  interior  of  the  sail  or  the 
appendages,  are  marked  for  exclusion  for  the 
computation  and  are  termed  “hole  point”.  Points  that 
surround  the  hole  points  become  new  inter-grid 
boundary  points,  which  are  called  fringe  or  receptor 
points,  and  require  boundary  values  to  be  applied. 
The  boundary  values  required  by  a receptor  point  are 
provided  by  interpolating  from  a donor  grid  that 
overlaps  the  region.  The  definition  of  which  points 
are  hole  points,  receptor  points  along  with  their 
corresponding  donor  members  is  termed  domain 
connectivity  information  and  is  provided  by  an 
overset  assembly  code,  SUGGAR  (Structured, 
Unstructured,  Generalize  overset  Grid  AssembleR). 
SUGGAR  is  designed  by  R.  Noack  [Noack,  2005],  to 
produce  overset  composite  grid  assemblies  of 
structured,  unstructured,  and/or  generalized  grids  for 
node  centered  or  cell  centered  flow  solvers. 

THE  INTERFACE  LIBRARY  - DIRTLIB 

The  domain  connectivity  information  obtained  from 
SUGGAR  is  transformed  to  the  solver,  and  the  hole 
points,  the  receptor  points  are  properly  identified. 
The  solution  on  appropriate  locations  in  the  donor 
grids  is  extracted  and  multiplied  by  the  interpolation 
coefficient  to  produce  the  value  required  for  the 
receptor  locations.  The  interpolated  values  are 


applied  at  the  receptor  locations  as  Direchlet 
boundary  values.  The  interpolation  and  distribution 
processes  are  performed  by  the  DiRTlib  (the  Donor 
interpolation  Receptor  Transaction  library).  DiRTlib 
is  developed  by  R.  Noack[Noack,  2005],  The  library 
simplifies  the  addition  of  an  overset  capability  by 
reducing  the  amount  of  code  the  flow  solver 
developer  must  write. 

THE  EXPERIMENT  - FLOW  ABOUT  ONR- 
BODY-1 

The  submarine  model  used  in  the  experiment  is 
DTMB  #5484  with  a length  of  5.27  m,  a diameter  of 
0.47m.  The  center  of  gravity  is  located  at  2.35m  (five 
times  the  body  diameter)  aft  of  the  forward 
perpendicular  along  the  centerline  of  the  hull.  The 
sail  has  a NACA  0014  section  with  a chord  length  of 
0.57  m and  an  aspect  ratio  of  0.27.  The  leading  edge 
of  the  sail  is  located  at  x/L  =0.2.  The  appendages, 
two  rudders  and  two  stem  planes,  are  identical  in 
shape  and  size.  They  have  NACA  0018  sections,  with 
a chord  length  of  0.18  m and  aspect  ratio  of  1.2.  It  is 
configured  as  ONR-Body-1  [Fu,  Atsavapranee,  Hess, 
2002].  The  tow  speed  is  5 knots.  The  corresponding 
Reynolds  number,  based  on  the  body  length,  is  13.5 
million.  The  details  of  the  test  procedure  and  setup 
can  be  found  in  reference  [Atsavapranee,  Forline, 
Fry.,  Hamiliton,  Percival,  Sung,  2004] 

By  using  stereo  particle  image  velocimetry 
(SPIV)  flow  patterns  with  all  three  velocity- 
components  around  the  model  were  gathered  and 
analyzed  for  three  model  configurations:  bare  hull, 
bare  hull  with  sail,  and  fully  appended.  The  main 
purpose  of  the  tests  is  to  correlate  to  the 
measurements  of  flow-field,  pressure  distribution  and 
overall  forces  and  moments  on  the  model  to  explain 
the  effects  of  the  dominant  flow  features  such  as  the 
sail  tip  vortex  on  the  flow  development  and  the  cross 
flow  separation  around  the  hull.  Test  data  are  well 
documented  and  are  suitable  for  the  validation  of  the 
newly  modified  and  assembled  UNCLE-overset 
code.  At  present,  only  a portion  of  the  test  matrix  will 
be  numerically  simulated. 

GRID  CONSTRUCTION 

The  ONR-Body-1  contains  all  the  essential 
components  of  an  operational  submarine  except  the 
propulsor.  namely,  the  hull,  the  rudders  and  the  stem 
planes.  The  construction  of  the  overlapping  grids  for 
the  computational  domain  is  greatly  simplified  since 
there  is  no  requirement  to  match  the  interface  grid 
points,  only  sufficient  overlapping  between  grids  is 
necessary. 


Based  on  the  same  geometrical  database, 
each  grid  about  each  individual  component  can  be 
generated  independently.  The  Hull  surface  is  first 
discretized  with  a network  of  four  sided  surface 
panels.  The  nodes  of  the  panels  lie  on  the  surface 
definition.  The  panel  density  can  be  varied  to  suit  the 
local  surface  curvature  or  the  intersections  with  the 
other  components.  The  body-conforming  volume 
grids  are  then  created  from  the  surface  grids  by 
hyperbolic  marching  [Chan,  1993],  The  marching 
process  can  be  control  by  specifying  the  off  body 
distance  of  the  first  grid  point  and  the  expansion 
ratio.  The  hull  grid  is  subdivided  into  a set  of 
overlapping  grids,  for  parallel  processing.  A simpler 
and  courser  far  field  grid  in  the  spherical  form  that  is 
overlapping  with  the  hull  grid  is  generated  to  extend 
the  far  field  boundary  to  three  body  lengths  upstream 
and  downstream.  The  sail  grid  itself  is  made  of  a set 
of  overlapping  grids.  It  has:  (1)  a collar  grid  that 
covers  the  sail/hull  intersection  region,  (2)  a wing 
grid  that  covers  most  of  the  sail  region  and  (3)  a tip 
grid  that  wraps  aroimd  the  tip  of  the  sail.  Such  grid 
topology  was  suggested  for  wing/fuselage 
configurations  in  reference  [Roger,  Roth,  Cao, 
Slotnick,  Whitlock,  Nash,  Baker,  2000],  Once  the 
sail  grid  is  completed,  it  can  be  added  to  the  hull  grid 
to  form  the  sail/hull  configuration.  This  procedure 
can  be  repeated  by  adding  the  grids  for  mdders  and 
stern  planes  to  the  hull  grid  to  form  computational 
grids  for  a fully  appended  configuration. 

Upon  the  application  of  the  grid  assembling 
code  SUGGAR,  the  hull  grid  points  that  lie  within 
and  immediately  surrounding  the  sail/appendages  are 
removed,  and  the  fringe  points,  donor  points  and  the 
hole  points  are  identified.  Furthermore,  the 
interpolation  stencils  associated  with  each  fringe 
point  are  computed  and  stored  as  part  of  the 
connectivity  information. 

Figure  I shows  the  hilly  appended  ONR 
Body-1.  Notice  that  there  are  no  fillets  in  the  juncture 
regions. 


Figure  2 shows  the  overlapping  surface 
grids  at  the  sail/hull  juncture  region.  It  can  be  seen 
that  the  sail  grid  is  composed  of  a collar  grid,  a wing 
grid  and  a tip  grid.  Similarly,  Figure  3 shows  the 
overlapping  surface  grids  at  the  appendages/hull 
juncture  region.  Clearly  all  the  surface  grids  are  near 
orthogonal  and  well  define,  and  at  the  overlapping 
regions,  the  grids  of  different  components  are  about 
the  same  size. 


Figure  2:  Overlapping  surface  grids  at  sail/luill 
juncture  region. 
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Figure  3:  Overlapping  surface  grids  at 
appendages/hull  juncture  region. 


Figure  1:  Fully  appended  ONR  Body-1 . 


Figure  4:  Complete  overlapping  surface  grids  on 
ONR  Body- 1 . 


Figure  4 shows  the  complete  surface  grids 
for  the  frilly  appended  configuration.  The  grid  density 
is  higher  near  bow  and  stem  of  the  hull  so  that  the 
surface  curvature  is  well  represented.  For  the  same 
reason,  the  grids  are  clustered  near  the  leading  and 
trailing  edges  of  the  sail  and  appendages. 

Figure  5 shows  the  contour  of  the  local 
surface  panel  area  of  the  sail/hull  configuration,  the 
length  scale  is  normalized  by  the  hull  length.  The 
contour  reflects  the  distribution  of  the  surface  grid 
density;  surface  panel  area  is  smaller  at  the  region 
where  the  radius  of  curvature  is  small. 
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Figure  5:  Contour  of  local  surface  panel 
area. 

Figure  6 shows  that  the  sail  and  appendage 
grids  embed  in  the  hull  grid.  Figure  7 shows  that  the 
hull  grid  embeds  in  the  far  field  grid. 


Figure  6:  Sail  and  appendages  grids  imbedding  in 
the  hull  grid. 


Figure  7:  Hull  grid  embedding  in  far  field  grid. 

The  total  number  of  grid  points  is  about  10 
millions  for  a fully  appended  configuration.  For  the 
purpose  of  multi-processing,  the  grid  systems  are 
divided  into  blocks;  each  block  has  approximately 
quarter  million  nodes.  Twenty-two  and  thirty-four 
blocks  are  used  for  the  sail/hull  and  the  frilly 
appended  case,  respectively. 

NUMERICAL  SCHEME 

The  numerical  solver  used  in  present  simulations  is 
based  on  the  CFD  code  UNCLE  developed  at  MSU. 
It  is  based  on  a three-dimensional  incompressible 
Navier-Stokes  equation  expressed  in  the  conservative 
law  form.  The  equations  are  coupled  with  the 
continuity  equation  through  the  artificial 
compressibility.  The  diffusive  tenns  are  discretized 
with  a central  differencing  scheme,  and  the 


convective  terms  are  discretized  with  the  Roe’s  flux- 
difference  splitting  scheme  [Roe,  1981],  For  the 
steady  state  simulation,  a first  order  Euler  backward 
approximation  is  used  to  represent  the  time 
derivative.  To  accommodate  the  rotational  motion, 
the  equations  of  motion  are  first  formulated  in  a 
rotating  coordinate  system  attached  to  the  moving 
body,  but  are  then  recast  in  the  absolute-flow 
variables.  Such  a formulation  is  popular  in 
turbomachinery  applications.  It  is  believed  that  such 
formulation  allows  more  accurate  calculation  of  flux 
in  the  finite  volume  method. 

The  Navier-Stokes  equations  of  the  three- 
dimensional,  time  dependent,  incompressible  flow 
can  be  written  in  a curvilinear  coordinate  and  in  a 
strongly  conservative  law  form  as 

[y ] + (E  - Ev  )e  + (F  - Fv )„  + (G  - Gv )?  = 0 (D 


assumed  high  and  approaching  infinite.  The  shear 
fluxes  (EV,FV,GV)  are  then  relatively  small  in 
comparison  with  their  corresponding  inviscid 
fluxes (E,  F,G).  It  can  then  be  assumed  that  the 
waves  that  travel  across  a cell  face  are  normal  to  the 
face.  Subsequently  the  three-dimensional  system  of 
partial  differential  equations  is  spatially  split  into 
three  quasi  one-dimensional  conservative  laws 

t]+h«=°.  <4> 

V J Jt 

where  9 is  8,  T],  or  C,  and  H is  E,  F or  G.  Defining 

computational  cells  with  their  centroids  at  m (m=i,  j, 
or  k),  and  their  cell  interfaces  at  m±l/2,  an  implicit 
discrete  approximation  to  Equation  4 can  be  written 
as 


Where  Q = (p,  u,  v,  w) T 
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with  the  inviscid  fluxes 

E*  = (pu,  u2  + p,  uv,  uw)T 
F*  = ({3v,  uv,  v 2 + p,  vw)T 
G*  = (pw,  uw,  vw,  w2  + p)T 
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.UjJ 

m m 

AQ”  + AH"  =0, 


(5) 


(2)  where  x is  time  step  size,  AQ"  = Qn"  -Qn  and 

AmH  = [Hm,„  - Hm_]/2]/A  9 and  D - 5H/dQ 

the  flux  Jacobian.  For  the  construction  of  an 
approximate  Riemann  solver  for  the  initial  value 
problem  in  Equation  5,  the  flux  at  the  cell  interface 
m+1/2  that  separates  the  right  from  the  left  traveling 
wave,  can  be  define  as 

(3)  H“+1/2  = + (A  Hm+1/2)n,  (6) 

or 


and  shear  fluxes 

Ev,Fv,Gv 

/?  is  the  artificial  compressibility  which  connects  the 
continuity  equation  with  the  momentum  equation,  J 
is  the  transformation  Jacobian  which  relates  the 
physical  domain  to  the  computational  domain,  and 

8x > £y, > £z > fix > fly > flz > C x Xy > Cz are  transformation 
metrics. 

In  order  to  apply  Roe’s  concept  of  flux 
differencing  to  the  convective  and  pressure 
differentiation  terms,  the  Reynolds  number  is 


h:+1/2  = h:+1-(ah;+1/2)",  (7) 

or  the  simple  average  of  Equations  6 and  7.  At 
present  calculation.  Equation  6 is  used,  since  it  has 
less  operation  count.  In  Equations  6 and  7,  the 
superscript  - and  + of  the  flux  differencing  term 
represents  the  right  and  the  left  traveling  wave  across 
the  interface  m+'A,  respectively.  That  is 

(ah„„)"  = [ra-r-'(Q1,-Q^)F  (8) 

and 


(m;  J = [ra*r-'(qL  -q;  )]”.  (9) 


with  R,  R 1 represent  right  and  left  eignvector  of  the 
eigen  system  of  flux  Jacobian  D,  respectively, 

A ' , A represent  positive  and  negative  eigen  value 
of  the  system,  respectively,  and  Q^,|,Q„  represent 
the  dependent  variable  to  the  right  and  the  left  of  the 
interface,  respectively.  For  the  first  order 

approximation,  QV,  = Qm+1  and  Q|n  = Qm  ■ 

To  evaluate  Equation  6 numerically,  a Roe’s 
[Roe,  1981]  mean  value  matrix 

Dm+i/2=D(Qm.Qm+i)  is  introduced  with  the 
following  properties; 


^m+l/2^m+l/2Q  — AHm+1/2  5 


(10) 


and 


D(Qm,Qm+1)  - Dm,  if  Qm  — > Qm+i  • (11) 


The  conditions  in  Equations  1 0 and  11  are  to  ensure 
that  the  solution  to  Equation  5 is  conservative  in 
space.  For  the  incompressible  flow,  the  only  state  of 
the  dependent  variable  Q that  satisfies  these 

properties  is  the  simple  average  of  theQ^+1, 
and  Qt'n  . Accordingly,  the  first  term  on  the  right  hand 

side  of  Equation  6,  H "n , is  computed  with  the 
metrics  evaluated  at  the  interface  m+Vi  and  the 
dependent  variable  Q evaluated  at  centroid  m,  that  is 

Q=Qm  • f°r  ^c  second  term  on  the  right  hand  side 
of  Equation  6,  (A  Hm+1/2)n , the  eignvectors  and  the 

eignvalues  R,  R 1 , A , arc  evaluated  with  the 
metrics  and  the  dependent  variable  defined  at  the 
interface  m +!4,  that  is  Q=  j (q*+1  + Q!n). 

To  enhance  the  accuracy  of  the  solution. 
Van  Leer’s  [van  Leer,  1974]  MUSCL  (monotone 
upstream-centered  scheme  for  conservation  laws) 
approach  is  adopted.  The  dependent  variables  to  the 
right  and  to  the  left  of  the  interface  of  the  cell  are 
reconstructed  as 
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R 

m+1/2 


Q 


m+l 


iL 

4 


_-(l  + KXQm+.  -Q„)_ 


and 
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L 

m+1/2 


0 -KXQm  — Qm  l)  ' 
_+(i  + kXQ„„i  -qJ_ 


(12) 


(13) 


If  cp  =0,  the  Equations  12  and  13  recover  the 
first  order  approximation.  If  cp  =1  and  K=  -1,  the 
approximation  is  second  order.  If  cp  =1  and  K=T/3, 
the  approximation  is  third  order.  An  alternative 
approach  to  the  reconstruction  is  to  apply  the  stencil 
shown  in  Equations  12  and  13  to  the  characteristic 
variables  W in  the  computational  space, 

whereW  = R 'Q . Once  Wn^+1/2  and  W^+1/2 arc 
obtained,  the  independent  variables  Q^+1/2  and 
Qmii/2  can  be  recovered  by  pre-multipling 

W,n+1/2  and  Wi,/,  with  the  right  eignvector  R 

computed  at  the  interface  m+1/2.  Tliis  reconstruction 
process,  identified  here  as  weighted  MUSCL  here, 
requires  higher  operation  count,  but  it  takes  the 
coordinate  non-uniformity,  stretching  and  orientation 
into  account,  since  the  characteristic  variables  are 
function  of  grid  position  and  matrix.  The  scheme 
described  here,  has  been  used  to  compute  the  flow  in 
linear  cascade  with  tip  clearance  [Yang,  2002], 
Similarly,  a fifth  order  stencil,  suggested  by  Rai  [Rai, 
1987]  and  used  by  Rogers  and  Kwak  [Rogers,  Kwak, 
1990],  is  derived  and  applied  to  the  characteristic 
variables  here.  That  is 


W*+1/2  = wm+] 

+ a)  [+2Am+5/2W  - 1 1 Am+3/2W 
■24Am+1/2W  + 3Am_1/2W] 


Am-l/2 


and 


(14) 


+ i[-2Am_3/2W  + llAm_1/2W 
+ 24Am+1/2W  - 3Am+3/2  W] 


(15) 


The  weighted  MUSCL  schemes  described 
here  are  implemented  in  the  UNCLE  code  for  present 
calculations. 

The  numerical  viscous  fluxes 
defined  in  Equation  (1  can  be  obtained  through 
central  differencing  technique  without  using  the  thin- 
layer  approximation.  The  differencing  equations  that 
represent  the  frill  Navier-Stokes  equation  formulated 
in  conservation  form  can  then  be  assembled.  The 
system  of  equations  is  solved  with  Newton  iteration 
method,  q-to  type  turbulence 

NUMERICAL  RESULTS 

Since  the  grid  systems  on  solid  surface  are 
overlapping,  direct  integrating  forces  over  the  surface 
would  result  in  the  overlapping  regions  being  counted 
more  than  once,  special  attentions  are  needed.  At 
present,  the  force  integration  routine  FOMOCO 
[Chan,  Buning,  1996]  was  used.  The  FOMOCO 
utility  is  a stand  alone  software  package  for 
computing  flow  coefficients  (force,  moment,  and 
mass  flow  rate)  on  a collection  of  overset  surfaces 
with  accurate  accounting  of  the  overlapped  zones. 
Tire  computational  results  presented  here  are  for  tire 
cases  of  flow  over  the  hull/sail  combination  and  fully 
appended  configuration  at  various  yaw  angles.  To 
facilitate  the  study  and  the  discussion  of  the  forces 
acting  on  the  body,  the  body  is  divided  into  several 
zones,  those  are,  the  bow,  the  nrid-body,  the  stem,  the 
sail  and  the  appendages  as  shown  in  Figure  8.  The 
bow  section  extends  from  the  forward  perpendicular 
to  x=0.129L,  from  there,  the  nrid-body  extends  to 
x=0.618L.  The  section  beyond  x=0.618L  is  defined 
as  stem  section. 


Figure  8:  Domain  division  for  force  and  moment 
integration. 


At  the  zero  degree  angle  of  attack,  the  axial 
velocity  contour  along  the  hull  for  the  sail/hull  and 
the  fully  appended  configurations  are  shown  in 
Figure  9.  The  axial  locations,  where  the  axial  velocity 
contours  are  shown,  are  x = 0.4L,  0.5L,  0.6L,  0.7L, 
0.8L,  0.9L  and  0.99L,  respectively.  Figure  9 shows 
that,  at  the  straight-ahead  condition,  the  axial  velocity 
profile  upstream  of  x=0.9L  is  not  affected  by  the 
present  of  appendages.  The  shape  of  the  contours 
suggests  the  existence  of  a horseshoe  vortex  behind 
the  sail.  Tire  turbulence  coordinate  y+  of  the  first  off 
body  grid  point  over  the  body  surface,  for  the  frilly 
appended  case  at  Reynolds  number  at  13.5  million,  is 
shown  in  Figure  10.  Over  the  majority  of  the  surface, 
the  value  of  y+  is  under  0.005.  This  is  an  indication 
that  there  is  a sufficient  number  of  grid  points  near 
the  body  surface  for  the  q-co  turbulence  model  to 
resolve  the  sub  layer.  Forces  integration  over  the  bow 
sections  for  the  sail/hull  and  fully  appended 
configurations  produces  identical  results.  The  same 
conclusion  is  observed  for  the  mid-body  sections  and 
the  sail  sections.  In  comparison,  the  stem  section  for 
the  fully  appended  configuration  produces  slightly 
more  drag,  because  it  has  less  pressure  recovery. 
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Figure  9:  Axial  velocity  contour  along  the  hulls  at  0° 
angle  of  attack. 
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Figure  10:  The  turbulence  coordinate  y+  for  the  fully 
appended  configuration  at  0°  angle  of  attack. 


Figure  11:  Law  of  wall  plot  at  hull,  X/L=0.365 

The  axial  velocity  profile,  in  term  of  the 
turbulence  coordinates  y+  and  u+,  at  the  mid-body 
section  of  the  hull  where  X/L  = 0.365,  is  shown  in 
Figure  11.  At  Reynolds  number  13.5  million,  there 
are  about  10  grid  points  within  y+=1.0.  The  present 
grid  system  is  suitable  for  computations  at  a much 
higher  Reynolds  number.  The  constant  4.9  associated 
with  the  logarithmic  function  in  Figure  1 1 is  for  the 
smooth  surface. 

The  computed  total  drag  coefficient  for  the 
sail/hull  configuration  is  0.000800.  The  contributions 
from  the  bow,  the  mid-body,  the  stem  and  the  sail  are 
14.375%,  51.00%,  29.50%,  and  5.125%  of  the  total, 
respectively.  The  pressure  drag  contributes  9.875% 
of  the  total  and  the  viscous  drag  contributes  90.125% 
of  the  total.  The  computed  total  drag  for  the  fully 
appended  configuration  is  0.000915.  The 
contributions  from  the  bow,  the  mid-body,  the  stem, 
the  sail  and  the  appendages  are  12,677%,  45.136%, 
28.087%,  4.481%  and  9.727%,  respectively.  The 


pressure  drag  contributes  13.989%  and  the  viscous 
drag  contributes  86.011%  of  the  total. 

Figure  12  shows  the  axial  velocity  contour 
along  the  fully  appended  hull  at  12°  yaw.  The  axial 
locations  for  the  contours  are  x/L=0.3,  0.4,  0.5,  0.6, 
0.7,  0.8  and  1.0.  The  contours  indicate  the  passage  of 
the  vortex  shedded  from  the  sail  tip,  the  primary  and 
the  secondary  separation  vortices  from  the  hull 
surface.  Figure  13  shows  where  the  streamlines  emit 
from  the  separation  lines  on  sail  and  hull  surfaces. 
The  sail  tip  vortex  induces  a hull  bound  circulation 
that  has  equal  strength  but  opposite  in  direction. 
When  this  circulation  interacts  with  the  cross  flow 
caused  by  the  yaw,  a downward  force  acting  on  the 
stem  is  generated.  This  induced  force  causes  the  bow 
to  pitch  upward.  When  the  yaw  angle  increases,  the 
flow  on  the  leeward  side  starts  to  separate.  There  are 
two  primary  separation  lines,  one  upper  and  one 
lower  as  shown  in  Figure  13.  Due  to  the  interaction 
with  the  hull  bound  circulation,  the  strength  of  upper 
separation  vortices  becomes  weaker  and  the  strength 
of  the  lower  separation  vortices  becomes  stronger. 
The  induced  circulation  of  the  lower  separation 
vortices  in  turn  weaken  the  hull  bound  circulation 
caused  by  the  sail  tip  vortex.  As  a result,  the  strength 
of  the  normal  force  acting  on  the  stem  section 
gradually  reduced.  As  the  yaw  angle  increases 
further,  the  resulting  nonnal  force  on  stem  changes 
sign  and  causes  the  bow  to  pitch  downward 

The  force  and  moment  results  are  presented 
in  a coordinate  system  whose  origin  is  located  on  the 
forward  perpendicular  on  the  body  axis.  The  x-axis  is 
positive  pointing  downstream.  Tire  y-axis  is  positive 
pointing  portside  and  the  z-axis  is  positive  pointing 
downward.  This  is  the  coordinate  system  used  in 
reporting  the  experimental  data.  Figure  14  to  Figure 
19  show  the  comparisons  of  the  six  components  of 
the  force  and  moment  coefficients.  The  comparisons 
are  very  good  except  for  the  force  and  moment  at  the 
out  of  plane  component  at  higher  yaw  angles.  This 
can  be  caused  by  the  presence  of  the  two  stmts  that 
are  supporting  the  hull  from  its  side.  The  stmts  can 
generate  a circulation  on  the  hull  that  interacts  with 
the  hull  bound  circulation  induced  by  the  sail  tip 
vortex.  The  interaction  may  change  the  normal  force 
on  the  hull  when  the  yaw  angle  is  sufficiently  high. 
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Figure  12:  Axial  velocity  contour  along  the  fully 
appended  hull  at  12°  yaw. 
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Figure  13:  Streamlines  emit  from  the  separation  lines 
on  sail  and  hull  surface. 
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Figure  15:  Comparison  of  lateral  force  coefficient. 


Figure  16:  Comparison  of  normal  force  coefficient. 


Figure  17:  Comparison  of  roll  moment  coefficient. 


Figure  14:  Comparison  of  longitudinal  force 
coefficient. 
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